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ABSTRACT: The interfacial structure formed in thermoreversible associating polymer solu- 
tions is studied within the density functional approach based on Flory's arguments of tree-like 
configurations of cluster associations. The unique characteristics of the interfacial region can 
be described in terms of the monomer density along the interface. For a certain value of the 
association parameter, which controls thermoreversible chemical reaction between associat- 
ing functional groups, the density profile is not smooth and undergoes a sudden jump at the 
point where the number of bonded functional groups is small. Analytical expression for the 
interfacial tension is given and results of numerical calculations are presented. 
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1. Introduction 



The systems containing molecules capable of forming thermoreversible chemical bonds with 
each other are of significant interest from both fundamental and practical points of view. 
These types of systems are often referred to as weak gels0 because the association between 
molecules under certain physical conditions, in particular lowering the temperature, leads 
to the formation of a thermoreversible infinitely large cluster (network or gel). At higher 
temperature, the molecules only form small, separated clusters of finite size, in the so-called 
"sol" solution. 

The thermodynamical properties of weak gels have been modelled previously with various 
theoretical approaches and approximations.il"0 An extensive theoretical study of the sol- 
gel transition phenomena was recently presented in Ref. ^ where many existing theoretical 
models and treatments were surveyed in details, showing that the pregel regime (or sol 
fraction) can be well accounted for with the use of a mean field approach. The original idea 
was first introduced by Floryi and Stockmayei0 who assumed that all finite-size cluster 
associations in the system consist of tree-like branches only. The analytical extension of 
this approach to the postgel regime, where gel fraction appears, gives the classical theory 
of gelation now widely accepted.0 The sol-gel transition within the frame work of this type 
of theory demonstrates smooth transition between the two states in the thermodynamic 
behavior, without any singularities typical in terms of phase transitions. 



However, as was first discussed in Ref. 12, it is important to include the possibility 



of formation of closed cycles in the model, in order to describe an infinitely large cluster 
properly. As the result of such inclusion, the sol-gel transition has the characteristics of a 
first-order thermodynamic phase transition.! 

Phase behavior of associating polymer solutions was the target of these studies and no 
attention has been paid to their interfacial properties. In the present paper we study the 
interfacial properties related to the region between the two coexisting phases of associating 
molecules. As a simple model system we take the system of /-functional monomers, each 
carrying / functional groups A that can react with each other forming reversible chemical 



bonds A — A (see Fig. |l|a). To obtain the free energy of this model system we use the 
classical density functional description which deals with tree-hke clusters only (see Fig. |I]b). 
The paper is organized as follows. In Section ^ we give an explicit expression for the free 
energy functional of inhomogeneous solution of tree-like associating molecules. In Section 
^ we describe the analytical and numerical procedures for minimization of the free energy 
functional and derive the expression for the interfacial tension. Readers who are not inter- 
ested in technical details might skip sections |^ and ^, and go to section ^ directly where we 
discuss the obtained results. 

2. Free energy of an inhomogeneous system within the functional density ap- 
proach 

Consider a solution of /-functional monomers distributed in space with a nonuniform 
monomer density p(r). The free energy of this system can be divided into two partJll 

F([p(r)],T) = F*([p(r)],T) + F,,.([p(r)],T), (1) 

where the first term on the right-hand side contains the contribution from the excluded- 
volume interaction between monomer units and the second term is related to the structural 
free energy of the system of an ideal solvent. The excluded- volume interaction, in particular, 
can be accounted for by the following expression, which was proposed by Flory0 based on 
a lattice model consideration, 

F*([p(r)], T) = ^ /{[I - vpir)] ln[l - ^;p(r)] - xv'p\r)}dr, (2) 

where T is the temperature of the system in ks units, v is the volume of one lattice cell (the 
volume of one monomer) and x is the Flory-Huggins constant, with xT being the energy 
gain per monomer- monomer contact. 

Our main purpose in this section is to introduce the structural term Fstr{[p{T^)],T) as 



derived within the density functional approach in Ref. 0,|T^. Assuming that the average 
distance tq between monomers is much less then the average length a of a chemical bond 
A- A. 



ro -C a, (3) 
we write the structural free energy in the form 

Fstr{[p{r)],T)/T = -^,4p(r)] +min{5eo™4p(r),PA(r)] - SA[pA{r)]}. (4) 
Here the first term 

SUpir)] = - [pir)\n(^-^]dr (5) 



/! 

corresponds to the translational entropy of the monomer units in space. The minimum of 
the second term is considered in the configurational space of all possible density distributions 
PA(r)- We also have 

ScomMr), PA{r)] = -J fp{r) [T{r) ln(r(r) + (1 - r(r)) ln(l - r(r))] dr, (6) 

where r(r) = p^(r)//p(r) is the local content of conversion, corresponding to the entropy 
of selecting reacted functional groups distributed with the density p^(r). 

The other term in eq. (^), Sa, corresponds to the entropy of thermoreversible chemical 
bonds. Here we assume that under condition (H) Sa depends only on the density of reacted 
functional groups and it does not depend on a particular monomer structure. Thus Sa can 
be calculated for the simplest system of fully reacted monofunctional units. From Appendix 
we have: 

pA{r){gip){r) 



where pyi(r) and '?/'(r) are dependent equation 



dr, (7) 



PA{r) = ^p{r){gip){r) = ^/j{r) / ^(r - r')^jj{r')dr 



with g{r) being the statistical weight of a bond with length r. For a homogeneous system 
we can rewrite eq. (|^ in the form 

SAiV, pa) =VP^ln^ (9) 
2 e 



where 

' g{r)dr. (10) 



The association caused contribution was also obtained by Semenov and Rubinstein in 
Ref. P who considered the homogeneous systems only. After minimization of the second 
term in eq. (^) we recover the expression which is widely used to consider phase behavior of 
associating polymer systems 

F^riT, V, p)/{TV) = pin (^] +fp(l+ ln(l - T)) , (11) 
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and 



fkp = j^,. (12) 

We shall not discuss details of phase behavior in this paper, but would like to mention that 
for certain values of parameters x and g it is more preferable for the system to be separated 
into two macroscopic sol phases with respectively different densities of monomers pi and p2 
which could be obtained from the equations 

Pi{T, pi) = p2{T, p2) (13) 
7ri(T,pi) = 7r2(T,p2) (14) 

where pi and tTj are the chemical potential and osmotic pressure of phase i {i = 1,2). We 
should also point out that if the density of monomers in the phase exceeds the critical value 

then the gel phase (infinitely large network) will be formed. 

In the next section we consider only the systems where two sol phases with different 
densities are present, describe the procedure that minimizes free energy functional ([l|), and 
derive the expression for the interfacial tension between the coexisting phases. 
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3. Minimization of the free energy functional and the interfacial tension 

In order to proceed further, we introduce a coordinate system so that the z = surface 
coincides with the flat interface separating the two sol phases. For a density variation whose 
typical correlation length is much greater than the basic Kuhn length associated with a single 
bond, the operator g can be approximated^ 

where the chemical equilibrium constant k is deflned in eq. (|TDp and a is the average length 
of a chemical bond. The free energy per unit surface can then be expressed as 

J'(<P{z),M^)AA{z)yz\ (17) 

where S is the area of the interface, Li and L2 are the linear dimensions of each phase normal 
to the interface, (p = vp is the volume fraction of monomers and we will refer to (pA = vpA as 
the volume fraction of reacted functional groups. The dot symbol represents the derivative 
with respect to z 

(IS) 




dz 

The density free energy JF, which is a functional of the local volume fractions of monomers 
and reacted functional groups, takes the form 

(j)A, </>a) = y^bis{(p) + ^comb{(p, (Pa) + ^a{(Pa, (Pa)- (19) 
The flrst term in ([l9|) corresponds to the free energy of the broken-link system 

J',is{(p) = 01n0 + (1 - 0) ln(l - 0) - X(P\ (20) 
the second term in (p!9| ) is simply combinatorial 

TcomM, 0a) = /0(r In r + (1 - r) in(i - r)) , r = (21) 

and the third one corresponds to the entropy of the chemical bonds 



9 A ^ S^A 

e 48 0A 



(22) 



where we have introduced a dimensionless chemical reaction constant 



g = k/v. 



(23) 



While the free energy across the interface is given by eq. ul7\ ), the surface tension is 
defined through considering the excess free energy in reference to the bulk phases. We can 
then write 

L2 



T 

a = — min 

V 



where the terms 



J (^(0, <j)A, <Pa) - -^1'°'") dz + J (^(0, <f)A, ^a) - J't"') dz 

-Li 



(24) 



^hom _ ^homi 



F. 



1 - r. 



(^ = 1,2), 



(25) 



are simply free energies of two homogenous phases with monomer volume fractions 0i and 
02 respectively on the left and right sides of the interface. 

The minimization of the functional (0) is considered with the constraint 



L2 
{(j){z) - (j)i)dz + j (0(z) - (l)2)dz = 

-Li 

and the volume fractions 0i and 02 are obtained from eqs. ([ 
explicitly, 



(26) 



(|T^, or, written more 



horn 



hom 



hom 



^hom ^hom 



(27) 
(28) 



Since any contributions to the free energy density which are independent of 0, or linear 
in 0, could be dropped or added,0 we can introduce a new free energy density function 

^(0,0A,0a), 



^(0, 0^) = ^(0, 0^, 0a) + A0 + S (29) 

where the parameters A and i? are to be determined from equations 

+ A(j)i + B = (30) 
^•2''°'" + A(j)2 + B = (31) 

We can now rewrite eq. (E^ in the form 



+00 



a = — min ^ / ^"(0, 0^, 0^)^^ ^ (32) 



V 

where we have rewritten the integration hmits in terms of infinity since JF approaches zero 
far from the interface. 

According to a standard minimization procedure we introduce two Euler's equations 

^-^^ = (34) 

d(pA dz d(t)A 

As the function T is independent of a first integral can be used to replace eq. (|3 



^-0A^ = O. (35) 



From eq. (|33D 



+ v4 + /ln(l -r) = (36) 

UlfJ 

one can find the equilibrium profile of conversion F 

r„).l-exp(-i^?|.-^). (37) 

This expression is valid across the entire interface and it gives a relation between conversion 
F and monomer volume fraction 0. 
From eq. (^) it follows that 

= + A0 + 5 + /0(F In F + (1 - F) ln(l - F)) - ^ In ^ (38) 
48 0A 2 e 



Taking into account eq. (^) and the definition of F (see eq. (pT])) we can now rewrite eq. 
(0) for the interfacial tension 

2T T 

CT= — J J^„{(t){z))dz, (39) 



-oo 



where 

:FM = ^Us{<P) + A<p + B + f<p In ^^^(fir(^)p + ln(l - r(0))) (40) 

Furthermore, 0(2;) can be written as 

= /(0F + 0f ) = /0 (^F + <j)^^ (41) 

Finally, with the help of eqs. (^) and (^OD, 0(2;) obeys a simple differential relation 

•2 48 0FJ-, 



a2 



(42) 



where 



aF 1 a^JS 



j-^(i-r). (43) 



Therefore, the volume fraction profile (j){z) can be obtained using the simple iterative 
numerical procedure 

0(z,+i) = 0(z,) + ^{(f){zi)){zi+i - Zi) (44) 



where (f){(f)) is defined by eq. (^2|). The interfacial tension a can then be easily found through 
eq. (|39|) after the calculation of the density. 



4. Results and discussions 

In this paper we study the temperature dependence of the interfacial tension in the solution of 
associating monomers capable of forming thermoreversible chemical bonds with each other. 
The functionality / of a monomer unit defines the structure of associations that could be 
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formed . If / = 2 then the system would contain hnear molecules of various lengths. For 
/ = 3 more complex structures could be formed as shown in Figure |l|. 

The temperature dependence of the Flory parameter x is considered to have the simple 
form 

Q , ^ 

X = (45) 

although more exotic temperature dependence could be introduced. We also introduce pa- 
rameter g defined through eq. (|23| ) in the form 

5- = fl'oexp(-Ex), (46) 

where is the chemical reaction constant at high temperatures and E9/2 is the energy of 
a chemical bond formed between two functional groups. 

As the first step, we consider the interfacial tension in a system of bi-functional units 
(/ = 2). If the chemical reaction constant g is high enough then all functional groups would 
be reacted and the system would contain one infinitely large linear chain. In this limiting 
case our theory should recover the treatment for the interfacial tension proposed by Helfand 
and TagamillUll. Writing free energy of the solution containing infinitely long linear chain 
in the form0 



V 



(V0)^ + (1 - 0) ln(l - 0) + x0(l - 0) 



dv (47) 



24 

one could obtain the expression for the interfacial tension ctq corresponding to the case where 
the rich phase (containing all chain monomers) and poor phase (only solvent molecules) are 
independently formed 



ao = — (^^ln(l-0) + x(l-</))-/^oJ (48) 

Here TyUg is the chemical potential and 0o is the volume fraction of monomers in the rich 
phase. These two parameters could be found from the equations 

/io = x(l-20o)-ln(l-0o)-l (49) 
/^o0o = X0o(l - 0o) + (1 - 0o) ln(l - 0o) (50) 
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which balance chemical potentials and osmotic pressures of bi-functional monomer units in 
the poor and rich phases. 

Figure H shows the temperature dependence of interfacial tension ao given by eq. ( ^Sf ) 
(open circles) together with a based on eq. (|39D for E = and different values of interaction 
parameters go introduced in eq. (|^). For large go (In go = 20, dotted line) we see that 
our result agrees exactly with the Helfand-Tagami theory thus confirming the model in the 
(/ = 2, (7 > 1) limit. 

Our model also accounts for the situation where a polydispersed system of linear chains 
system can be formed for smaller go- For smaller values of go [In go = 2, dashed line, 
and In go = 0, solid line) when the reaction of functional groups is incomplete one can see 
the singularities that might occur in the o"(T) behavior. Distinct singularity points can be 
identified where sudden changes in curvatures occur. 

Figure |^ demonstrates the interfacial tension of the / = 3 case where the branched 
associates may coexist across the interface and one can also observe the similar singularity 
behavior for cr(T). 

The described types of singularities in the interfacial tension behavior arise due to the 
singularities that can occur in the behavior of the density profiles (j){z) and </>a(-2) along the 
interface (smooth profiles without any singularities are shown at the top of Fig. |^). For a 
certain value of parameters x ^''^^ 9 o^i^ can expect a sudden change in the slope of (f)A{z) 
as shown in Fig. ^ This picture is observed if T + (p dT / d(f) introduced in eq. (^21) becomes 
negative in the interval (0',0") as shown at the bottom of Fig. |^. According to eq. (^) the 
slope of (j){z) becomes infinite at z' and z" which correspond to 0' and 0" and the slope of 
0yi(z) is negative in the entire interval [z', z"). 

The more complex behavior is observed when the equilibrium conversion r(0) defined 
by eq. (|37| ) becomes negative, as shown at the bottom of Fig. |^ for the interval (0",0"'). 
Physically F is a positive quantity and therefore the region (0", (f)'") should be skipped when 



calculating a defined by (|39D. Skipping this {(f)", (j)'") interval leads to the sudden jump of 
the volume fraction profile (j){z) at the point z" as shown at the top of Fig. |^. Similar 
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noncontinuous behavior of polymer density was previously discussed in Ref. ^ where the 
authors considered a polymer globule formed by a long polymer chain with saturated physical 
bonds. This type of behavior could be explained as follows. As one goes across the interface 
(see the top of Fig. ^) the volume fraction (p increases until it reaches the value 0" at the 
point z" . At the same time the number of reacted functional groups (or reduced density of 
reacted functional groups pa) decreases and becomes infinitely small as cf) reaches (f)" (at this 
point conversion F is zero). The fact that 0a = ensures that there would be no loss of 
conformational entropy for the system associated with the jump of the volume fraction at 
the point z" because square gradient term (see eq. (|3)) depends only on derivative of 0a 
which is continuous at this point. 

Hence, the behavior of interfacial profile cr(T) shown in Figure ^ could be explained 
as follows. As the temperatures of the system decreases and reaches Ti the interface first 
appears and the interfacial tension starts to increase. As temperature reaches the value of T2 
the change in the curvature of cr(T) profile is observed. This curvature change is connected 
to the appearance of the negative slope in the behavior of 0a(-2) (as shown in Fig. When 
the temperature of the system reaches T3 one observes the change in the slope of o"(T) profile 
due to the discontinuity of the monomer volume fraction 0(z) (see Fig. |^) which first appears 
at this temperature. 



5. Conclusion 

In this paper, we analyze the interfacial structure and calculate the temperature dependence 
of the interfacial tension a in associating solutions of /-functional monomers. We find that 
there are two types of singularities which occur in the behavior of the density profiles along 
the interface. At certain temperatures the density of reacted functional groups may decrease 
within some region inside the interface while the polymer density always increases. Also, 
the polymer density may exhibit a sudden jump at the point where the density of reacted 
functional groups becomes infinitely small. The first of these singularities leads to the change 
in the curvature of the cr(T) profile and the second one changes the slope of that profile. 
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A. Calculation of the entropy of thermoreversible chemical bond^iUli 

In this Appendix we cosider the system of fully reacted monofunctional monomers units. If 
p^(r) is the monomer density, then the entropy of chemical bonds Sa is defined as follows 



where Spair is the entropy of dimers formed due to the associations and Sid defined by eq. 
(^) corresponding to the translational entropy of monomers. We assume that external field 
(f{r) produces an equilibrium density PA(r); and in the mean field approximation 



Here -E[<y9(r)] and Fpair[f{r)] are the energy and free energy of dimers in the field ip{r) 




(Al) 



TSpair[pA{r)] = E[Lp{r)] - Fpair[ip{r)] 



(A2) 



respectively. They satisfy the following expressions! 








(A3) 




(A4) 



where N is the total number of dimers 




(A5) 



and 




(A6) 



(A7) 
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with g{r) corresponding to the statistical weight of a bond with length r. According to 
general thermodynamic relations^ the density of functional groups can be written 



as 



^^^^^ ^ SF,,,r[f{r)] ^ 2Ar>W(^^)(r) 



^(r)(#)(r) 



(A8) 



Using eqs. ( [A5[ ) and ( |A8| ) one can rewrite eq. ( |A4| ) in the form 



F,„„,Kr)] = -r/^ln 



Pa{t) f2 



dr 



(A9) 



Eliminating (f{r) from eq. ( [A3| ) with the help of eq. ( [A6| ) and substituting the result into 
eq. (|A2|) together with expression ( |A9|) we obtain 



5'pajr[PA(r)] 



In 



V'(r)p^(r) w2 



(ir, 



(AlO) 



and for the entropy of chemical bonds given by eq. ([Al|) we finally get 



PA{r) 



In 



PA(r)(#)(r) 
e-?/'(r) 



(All) 



which is the formula used in eq. (|^). 
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Figure captions 

Figure 1. Formation of a reversible chemical bond A — A between two three-functional 
monomers (a) and an example of a tree-like cluster formed by three-functional monomers 
(b). 

Figure 2. Interfacial tension a given by eq. (^) (t; = 1, a = 1) as a function of temperature 
for bi-functional units with f = 2, E = 0, and different values of Ingo. Open circles 
correspond to the limited case given by eq. (^). 

Figure 3. Interfacial tension a given by eq. (^) (f = 1, a = 1) as a function of temperature 
for three-functional units with / = 3 and values of interaction parameters shown in the upper 
right corner. 

Figure 4. Equilibrium profiles of volume fraction (solid line) and volume fraction of 
reacted functional groups (pA (dash line) along the interface in the case of / = 3 and values of 
interaction parameters shown in the upper right corner (top). Also shown are the equilibrium 
conversion F (solid line) and T + (j)dr/d(f) (dash line) as functions of volume fraction 0. 

Figure 5. See caption for Fig. |^. 

Figure 6. See caption for Fig. ||. 



17 



Figure A.Ermoshkin, I.Erukhimovich, J.Chen 




Temperature 



Figure A.Ermoshkin, I.Erukhimovich, J.Chen 




Temperature 



Figure A.Ermoshkin, I.Erukhimovich, J.Chen 




0.0- 

1 1 1 1 1 

-15 -10 -5 5 10 z/a 

Distance along the interface 




I I I I 

0.0 0.2 0.4 0.6 0.8 



Volume fraction 



Figure A.Ermoshkin, I.Erukhimovich, J.Chen 





0.0 0.2 0.4 0.6 0.8 

Volume fraction 



Figure A.Ermoshkin, I.Erukhimovich, J.Chen 





Volume fraction 



